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Abstract 

The integrable perturbation of the degenerate boundary condition (d) by the ipi^ boundary 
field generates a renormalization group flow down to the superposition of Cardy boundary states 
(+)&(—). Exact Thermodynamic Bethc Ansatz (TBA) equations for all the excited states are 
derived here extending the results of £Q to this case. As an intermediate step, the non-Cardy 
. . . boundary conformal sector (+)&(—) is also described as the scaling limit of an A4 lattice model 

, with appropriate integrable boundary conditions and produces the first example of superposition 

' of unitized Virasoro characters. 

(N 

' 1 Introduction 

On 



Two dimensional quantum field theories in presence of a boundary have received an increasing at- 
tention in the last few years because of their important role in condensed matter physics (especially 
CN ■ on quantum impurity problems) and in string theory (D-branes). Their integrable or conformal 
properties lead to the computation of a number of exact results (boundary S-matrices |2j, conformal 
Q ' boundary states [Sj, etc. See also jl] and references indicated there). An interesting problem is to 
connect different boundary conformal field theories (BCFT) by integrable renormalization group 
flows generated by boundary operators. 
£f) ' The tricritical Ising model is the conformal unitary minimal model with central charge c = 

7/10. Its phase diagram in presence of boundaries, with the bulk of the system maintained at the 
critical point, contains various fixed points (that are precisely boundary conformal field theories) 
connected by renormalization group flows jSlEl- Among them, those generated by perturbing a 
fixed point with the magnetic boundary field ip\^ are integrable. One of the two flows departing 
■ from the BCFT (d) moves to a superposition of Cardy-type boundary states, (+)&(—). This is the 
simplest occurrence of the general observation [7j that the boundary interactions do not make a 
distinction between pure states or superposition of states. Actually, increasing the central charge 
of the minimal models, one expect more and more the appearance of superposition of states |Sj so 
that the analysis presented here is a crucial step toward a quite general phenomenon. Physically 
speaking, boundary correlation functions of pure states satisfy the cluster property while this is 
not true for superposition of states. This is an indication of a first order phase transition occurring 
close to a BCFT that is not a pure state jH]- 

The flow (d) = (2,2) 1— ► (+)&(—) = (1,1) © (3,1) is studied here using a lattice approach 
to formulate TBA equations, as described with great detail in on which the present article is 
heavily built. For convenience, the same notations will be used. This flow was identified in |5j by 
the boundary S matrix; it preserves supersymmetry [2]; a TBA system for the ground state only 
was given jlUj in a different channel (periodic). 

The plan of the paper is as follows. In Section[2j the pattern of zeros for the boundary critical 
point (1, 1) © (3, 1) on a finite lattice will be described and the unitized partition function will 
appear as a superposition of unitized characters JJ. In Sectional the flow will be recognized to 
be in the "variable r" family, and the mapping between unitized characters will be provided. In 
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Section |1J the specific set of TBA equations for this flow will be derived and numerically solved. 
In Appendix A there is a summary of properties of Gaussian polynomials. In Appendix B the 
expression for the energy at the BCFT (1, 1) © (3, 1) will be proved. In Appendix C the selection 
rules for the integration constants will be given, that also apply to all the cases discussed in [Q. 

2 The boundary critical point (1, 1) © (3, 1) 

An N faces double row transfer matrix with boundary conditions (r^, Oj), i = 1, 2 on the left /right 
boundaries is indicated by a , Q2 (tt, £1,^2) where < u < A is the spectral parameter, A = ^ 
is the crossing parameter and £j are parameters associated to the boundary weights. If £ is chosen 
appropriately, the lattice boundary weight (r, a) scales to the conformal boundary condition [T^IT^] 
also labeled by (r, a) [II] . In the rest of the paper, the only case that will be used is: 

D>,6)=Dj 1|2jl (u,6,^A). (2.1) 

The parameter £2 will be kept fixed at the indicated value corresponding to the conformal boundary 
condition (2, 1). 

From [T3], it is known that the continuum scaling limit of the double row transfer matrix 
D iV (^, |A) is a realization of the (2,1) © (2,1) = (1,1) © (3,1) conformal boundary conditions. 
This is obtained if the boundary parameter is chosen inside the real interval £1 E [A/2, 5A/2]. The 
choice of the middle point of the interval will be useful later when moving off criticality by the 
introduction of an imaginary part for £1. The number of faces N is required to be even by the 
adjacency conditions. 

The pattern of the zeros of the eigenvalues D{u) of D(tt) is relevant to the classification of the 
states, both for the lattice and for the continuum theory, and must be determined with numerical 
observations performed on lattices of small size 1 N = 12 (in the name "numerics on D" was 
used) . As typical for the A4 lattice model , there are two analyticity strips 

A 3A 

strip 1: - - < Re(u) < — , strip 2: 2A < Re(u) < 4A 

and two types of strings. The 1-strings are single zeros appearing in the center of each strip 

Re(it) = — or 3A 

while the 2-strings are pairs of zeros (u,u') appearing on the edges of a strip, with the same 
imaginary part 

f (-A/2,3A/2), strip 1, 

(Re(«),Re(«')) = { 

{ (2A,4A), strip 2. 

In the upper half plane, let mi, W12 be the number of 1-strings and m, n2 the number of 2-strings 
in each strip. The following relations hold between the numbers of zeros: 

n x = to_ mi ^ ) 

> =>mi,m2 even. (2.2) 
n 2 = f -m 2 +1^0 J 

The states on the lattice and at the scaling limit are uniquely characterized by the non-negative 
quantum numbers {ij?^}; they have the topological meaning that for a given 1-string, 7^ i s the 
number of 2-strings with larger coordinate. This remains true at the scaling limit. 

lr The number of faces is mainly limited by the huge amount of RAM required to write the transfer matrix. 
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The determination of the scaling energy corresponding to each eigenvalue is the crucial element 
to distinguish states in the (1,1) or (3,1) sectors. Remembering that A^x = 0, A^i = ~, the 
separation between levels is quite sharp. In spite of this, the finite size effects for N = 12 faces are 
still dominant on the accuracy of the energies, apart for the first few levels E ^ 5, so the complete 
distinction of states among the two sectors requires informations from the TBA equations that will 
be introduced later. Anticipating here these results, one has: 

(1, 1) : lf ] < n 2 or m 2 = 0, (2.3) 
(3,1): lf ] =n 2 . (2.4) 

This means that the discrimination is given by the lower object in the second strip: if it is a 2-string, 
the sector (1, 1) is selected, if it is a 1-string the sector (3, 1) is selected. The given conditions imply 
two useful necessary conditions for a state to be in the corresponding sector: 

(1, 1) ^ lf ] < n 2 ^ - m 2 0, (2.5) 
(3,1) m 2 ^2^mi^2. (2.6) 

The former is an enforcement of the second expression in (|2.2|) . The energy corresponding to each 
state is finally given by 



i=l,2 k=l 



24 4 2 H >2 



C = ( _ „ I , 5 a h : Kronecker's delta, m = ( mi 
\-l 2 J ' \m 2 

C is precisely the A 2 Cartan matrix. The complete and independent proof of this expression will 
be given in Appendix B. 

Here a crucial observation can be done that there is a violation of the topological rule for 
the energy. Indeed, when only one sector is involved, if the relative position of a 1-string and a 
neighboring underlying 2-string are exchanged the energy increases of 1 unit. The opposite exchange 
reduces the energy of 1 unit [121 ■ Presently two sectors are considered and this rule is violated in 

precisely one case: if the lower objects in the second strip are a 1-string and an underlying 2-string 

t 2 \ 

(i.e. I\ = n 2 — 1), their exchange increases the energy of ^ units. The opposite exchange reduces 
the energy of ^ units. This corresponds to the Kronecker's delta term in (|2.7[) . In all the other cases 
the usual rule holds. Notice that this violation is expected from counting arguments, namely with 
the usual rule only, there is no way to match the states given by D(u) with the states contained in 
(1, 1) © (3, 1). In Tabled there is the first appearance of this violation, namely the infrared state 2 
(0 0|0 0), whose energy is 6, has two excitations, (0 0| 1 0) , with energy | + 5, and (1000|00), 
with energy 7. 

The finitized partition function is defined as a sum on the states of a finite lattice where the 
energy of each state is fictitiously assumed to be equal to its scaling limit energy: 

Z N (q) = q 24 ^ q * = q 24 ^ q h •*» (2.8) 

StatCS mi ,m 2 ,{/W} 

where q = exp(— 2irM/N) is the modular parameter. Now, it will be shown that the finitized 
partition function matches the sum of characters. The sum on states 1)2.8(1 can be split in two 

2 The notation indicates the content of zeros: for example, (3 2 2 1 0) is a state with 4 zeros in the first strip and 
2 in the second; their quantum numbers are respectively /J™ 1 = 3, I™ 1 = 2, . . .. If a frozen zero occurs 1121 . a parity 
a = ±1 is added. 
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parts, for l\ < n 2 (sector (1, 1)) and for I\ A> = n 2 (sector (3, 1)). Using the Gaussian polynomials 
and their properties as in Appendix A leads to: 

£ 



r(2) 



if l[ 2) < n 2 : 



r (2), r(2) 

q 1 ! T,,, '™2 



0<4^sS...</f ) <n 2 



m 2 + n 2 - 1 
m 2 



if /J 2) = n 2 



E 



r (2) , ,(2) 
fl + - 1 m 2 — qTl 2 



o<x^...crf ) =n 2 



m 2 + n 2 - 1 
m 2 - 1 



(2.9) 
(2.10) 



The sum on the first strip is not modified by this counting. The partition function now reads 



q 24 ^ q 

(1,1) 



\m T Cm 

4 



iri\ + n\ 




m 2 + n 2 - 1 


nil 




m 2 



(2.11) 



+ 24 



(3,1) 



rm T Cm— i ri2 



mi + ni 




m 2 + n 2 - 1 


mi 




m 2 — 1 



where the labels on the sums indicate that the sums on mi,m 2 are restricted by the constraints 
imposed by the corresponding sector (|2.2II2,6|) . The first sum is easily recognized to be the unitized 
character Xi^i{q)- I n the second term, the redefinition of m 2 — > m' 2 = m 2 — 1^1 and odd, leads 
to 

h mx m2 )c( m ^ + n 2 -\ = \{ mi m' 2 )ch) 
4 \m 2 J 2 4 \ m 2/ 

so that it becomes the finitized character X31 (<?) • This proves that the unitized partition function 
is the sum of the corresponding finitized characters: 



ZN(q)=X { 1 N 1 \q)+xZ'(D 



(AO, 



(2.12) 



and, consequently, that the pattern of zeros (|2.2j) and the distinction of sectors (|2.31 12. 4[) are 
consistent with the known scaling properties of the given transfer matrix. 



3 The flow: 3 mechanisms 



t I 



t i 
i 



X IX 31 

A 



I 2% 31 

B 



I 2X 3X 

c 



Fig. 1. The three mechanisms A, B, C respectively that change string content during the flow xi,i + X3,i 
*-* X2.2 which is the reverse of the physical flow. These mechanisms are illustrated for the states: A, 
(00|00) i ^ (00|0)+; B, (1 1|00) » (00|1)_; C, ( ) ^ ( |0)_. 



Considering the transfer matrix D(«, £i a tt) with £i at t = |A+z£/5, one can compute the expected 
endpoints, 

3 _ ^ | X2,2 if £ -> ±00 

I Xi,i + X3,i if C = 



(3.1) 



4 



Tabic 1. Flow Xi,i + X3,i l— * Xi,i (reverse of the physical flow). The explicit mapping of states from IR 
to UV up to the IR level 7 are presented here. n IR , n uv are the excitation levels above the ground states, 
respectively A^i = and A2,2 = gjj- The IR states in the two sectors are easily recognized because in (3, 1) 
they are separated by A.3 t i = | from the true ground state. 
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as it was done in (Tj for various types of boundary. The choice of the real part Re(£i at t) = |A 
makes the double row transfer matrix real symmetric on the center of each strip, Re(u) = 4, 3 A 
so that the eigenvalues and the energies are real. Any different value for the real part would lead 
to complex energies and non-unitary flows. Summarizing, one has integrability, correct endpoints, 
unitarity, and uniqueness so the identification of this lattice flow with the renormalization flow 
(2,2) — > (1, 1) © (3, 1) is completely natural. 

From extensive "numerics on D" one can figure out the correspondence between infrared (IR) 
and ultraviolet (UV) states, Tabled and the three mechanisms A, B, C active to change the pattern 
of zeros. They are the same known for the "variable r" flows in [Q so the present flow is part of that 
family. For their complete description the reader is referred to while here there is a summary 
of the properties and in Fig. ^ there is an example of their action. 

The mechanisms A,B,C force the following changes for the various parameters: 

m 2 R i — ^ mi Y = m 2 R - 1, (7=1. (3.2) 
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where no label IR or UV is used for parameters that do not change during the flow. 

An important consistency check is to show that the three mechanisms respect the counting of 
states at the endpoints of the flow. From (|3.21 13.31 13.4(1 one sees that for each IR state there is 
precisely one applicable mechanism so that the IR counting of states is complete. Moreover, using 
the recoursive identities among Gaussian polynomials (|A.21 IA.3|) the IR partition function in (|2.11j) 
naturally splits into six terms precisely associated with the three mechanisms, one for each sector: 
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(3.5) 



The labels (1,1), (3, 1), A, B, C indicate that the corresponding sums on mi, m!> R are restricted by 
the constraints imposed by the respective sector and mechanism. Notice that in the (1, 1) sector 

r(2) IR 



the quantum numbers in the second strip are bounded by I\ 
r(2) IR _ „m 



< n 



IR 



1 while in the sector (3, 1) 



the value I[ 



n 



is fixed and doesn't enter the combinatorics that is equivalent to say that 
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mJ^ — 1 zeros must be considered. The two sum constrained by A can be understood using (|A.5|) 

(1) (2) ^ 

for both the (/-binomial factors and correspond to summing under the constraint 4' = J f R = 0. 

Similarly, the sum on B becomes apparent using ()A.4|) for the first strip (/-binomial and (|A.5|) 
for the second strip (/-binomial and likewise the sum on C uses (|A,4|) for the second strip only. 

Two steps are now involved: the first is the mapping of energies, the second is the counting of 
states. The IR energy expression is given by ()2.7(l . For the UV sector (2, 2) is given here |12j : 
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(3.6) 



^2,2 = the rules for the various parameters are given in Table [21 An IR energy level at the base 
of a tower of states with string content fixed by (mi, mJ> R ) maps to a UV energy level according to 
the energy expressions (|2.7( 13.6 1) that hold at the two conformal endpoints of the flow 
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Using (|3.2() to (|3.4|) . the (/-binomials appearing in the IR partition function Q3.5J) are rewritten in 
terms of UV parameters. Also taking into account the mapping of the energies (|3.7j) one obtains: 
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The two terms corresponding to each mechanism can be summed with (|A.3|) so that the distinction 
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among sectors disappear 
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and one is left with the UV partition function, that is the (2, 2) finitized character. This argument 
shows that the IR counting of states with the three mechanisms induces the correct UV counting 
of states. 



4 The flow: TBA equations 

The TBA equations are obtained with the procedure given in [Q: the double row transfer matrix 
(|2.1|) eigenvalues can be normalized such that they satisfy a functional equation. The analytic 
content of the normalized transfer matrix eigenvalues is used to solve by Fourier transform the 
functional equation. As last step, a scaling limit is performed to obtain the equations for a contin- 
uum theory (TBA equations). The normalization factor is highly non-trivial in the sense that it 
contains both the spectral and the boundary parameter. The latter appears by a factor g(u, £iatt) 
that is common to all the "variable r" cases and, with reference to the center lines of the two 
analyticity strips and at the scaling limit, is given by: 

gi(x,£) = tanh^, (4.1) 
g 2 (x,0 = 1. (4.2) 

This procedure leads to the same TBA equations, energy expression and quantization conditions 
known for the "variable r" flows: 

rax (1) _ 

ei(:r) = -loggi(x,g)-E lQ g( tanh ^ 2 X )~K*L 2 , (4.3) 

k=l 

m 2 (2) _ 

e 2 (x) = 4e-"-logg 2 (x,e)-E lo g( tanh ^ 2 X )~K*Li (4.4) 

k=l 



where the integration kernel is given by K(x) = 27r * h and Lj(x) = log |1 + Sj exp(— ej)\. The 
denotes the convolution, (/ * g)(x) = dy f(x — y)g(y). The energy expression is 
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Tabic 2. Classification of the patterns of 1- and 2-strings at the endpoints of the flow. The parity a = ±1 
occurs in (2,2) because there are frozen zeros. The parities S\,S2 = ±1 occur in the TBA equations as 
integration constants (see Appendix C). The expressions for n\ are only used on a finite lattice because in 
the scaling limit n\ ~ N/2 — > oo. The number of faces in a row, N, is even. 
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and the quantization conditions are expressed in terms of the counting functions ^j{x) = 
-i€j(x - if): 

My { k ] ) = -^ 2 (yi 1} - i|) = 4 1} 7r, 4 1} = 2(4 X) +mi - k) + 1 - m 2) (4.6) 

MV®) = " »|) = 4 2)?r ' 4 2) = 2 (4 2) + "2 - fc) + 1 - mi. (4.7) 

Note the inversion of the indices: €\ is for strip 2 and e 2 for strip 1. This set of equations gives the 
complete description of the scaling energy for all the states of the flow under examination. The 
correct choice of the parameters mj, n 2 is suggested by the lattice (Section and a quick numerical 
check on the ground state leads to the values of the integration constants Sj. In summary, starting 
close to the IR point (£ — > +oo) and moving along the flow one has: 

mechanism A and B, mechanism C before the collapse: 

mi, m 2 even, ra 2 = ^ — m 2 + 1 ^ 0, si = s 2 = 1, 
mechanism C after the collapse point: 

= mo + 2, n" 2 \ , = n~ 2 \ = 1 — mi , the other variables remain unchanged. 

z rrtj— 1 

A small surprise is given at the IR point itself: performing the limit £ — > +oo on the TBA equations 
leads to two different setups for the two sectors, namely the previous summary is correct for both 
the sectors along the flow and for the sector (1, 1) only at the IR point. The sector (3, 1) at the IR 
point requires an odd number of 1-strings in strip 2 as well as s 2 = —1, as shown in Table EJ This 
odd behaviour of the flow close to the IR point is new: in all the previous cases, the movement 
off the IR point was observed to be smooth and the content of zeros never came up with a change 
in the value of m 2 . The reason can be easily found in the discontinuity of the following limit (eq. 
(3.70) of P): 

lim ipi(x)= lim ilog&(x,£) = { Z \ ! < T^' ( 4 - 8 ) 
x~^— oo x — > — oo I u it c; — +00 

In the first case, the odd parity of n k ([4.7)1 implies that a second strip zero can be found at 

(2) 

— oo. Indeed, a clear numerical evidence and the following simple argument show that y\ = —oo 
precisely for all the states in the (3, 1) sector. Looking at the difference between the endpoints of 
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ipi along the flow (assuming ym[ < +00) one has 

lim tpxfx) - lim Vi(^) = 7r(l + mi) = ir • "odd" (4.9) 

x— >— 00 x— >+oo 

so that the maximum number of zeros (1- and 2-strings) that find room in that interval 3 is [ mi 2 +1 j + 
1 = ^ + 1 that coincides with the actual number of 1- and 2-strings in the strip, mi + ni- So, all 

the locations are full and when 1^ = ri2, equivalent to the (3, 1) sector, the lowest location must 

be occupied by a 1-string. Observe that the alternative assumption ym\ = +00 leads to the same 

(2) 

conclusion because m± is even. Summarizing, along the flow one has y\ = —00; in the IR limit 
£ — > +00 this zero must leave the complex plane because of the second case in (|4.8|) and this leads 
to the rules for the TBA in (3, 1) sector, as in Table 

A numerical solution of this set of TBA equations for the first few states was done and is given 
in Fig. El 



5 Conclusions 

The lattice approach to boundary renormalization flows exposed in and extended here appears 
quite effective in obtaining the scaling energy for all the excitations. It should be interesting to 
understand something more about the role played by the picture of zeros in the continuum theory 
itself. Indeed, this structure survives the scaling limit and is relevant in the description of states. 

Moreover, there are structures that are common to the ^13 perturbations |15j of tricritical 
Ising model. They are the types of zeros and number of strips, corresponding to the number of 
pseudoenergies, the way these pseudoenergies are coupled to form the TBA system and the kernel. 
It should be nice to have track of this from the continuum theory itself. 

It is also an interesting open question to extend the present lattice approach to all the boundary 
flows for the whole series A-D-E of minimal models. 
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Appendix A. Gaussian polynomials 

The finitized characters are know to be expressed in terms of (/-binomials. They are combinatorial 
expressions defined by: 



m + n 
m 



n Ii Im—l 
/! =0/2=0 I m =0 



(A.l) 



Two recoursive relations are important for this paper: 



m + n 
m 



m + n — 1 
m — 1 



m + n — 1 

77! 



(A.2) 



m + n 
m 



m + n — 1 
m 



+ q r 



m + n — 1 
m — 1 



(A-3) 



' [x] means the integer part of x. 



10 




Fig. 2. Scaling energies for the flow X2.2 l_ > Xi,i+X3,i- The list of states is given in Tabled The intermediate 
region of the mechanism C levels (shown dashed) are schematic and have not been obtained from the solution 
of the TBA equations. 
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The following expressions are useful to understand the contribution of each mechanism to the full 
character. Indeed, the restriction I m > leads to 



n Ii Im-1 
Il=l/ a =l /m=l 

whereas the restriction I m = leads to 



i r' v 

n—l 1 l J m-1 



4=04=0 /^=o 



m + n — 1 
rri 



(A.4) 



n ii -?m-2 

EE- E 

/ 1= 0/ 2 =0 /m-l=0 



m + n — 1 
m — 1 



(A.5) 



Appendix B. Proof of the expression for the energy at (1, 1) (3, 1) 

The energy expression (|2.7|) is proved here using the continuum energy and the TBA equations in 
the limit £ — ► +00. The method is similar to the computation done in J2] at each critical point 
but, in the present case, one has to consider extra terms that appear because the limit £ — ► +00 
will be taken only at the end of the computation. 

_ (i) 

The first step is to transform 1)4.5(1 in a more explicit form. Its first term, e y k , can be 

expressed by the quantization condition 1)4.6(1 and then one adds = Y^k=i( n k % ~ ^(Vk )) ^ Iom 
(|4.7|) . Some numeric terms are obtained, given by sums and products of m« and integer numbers. 
They can be easily worked out with the help of the Cartan matrix defined in (|2,7[) leading to 



i=l k=l k=l 



+ 



"°° dx 
4tt2 



mi r / \ m 2 



EL\{x) v-^ L 2 {x) _ x 



The contribution from the integration can be worked out with the following procedure. The main 
point is to compute, in two different ways, the real integral: 

/+00 
dx [-e'jLj + Re(ej)4] (B.2) 
-00 

/+00 
dx [(log Sjtj)' log |1 + tj\ - Re(log(s/ i ))(log |1 + tj])') 
-00 

where tj = Sj exp(— ey). On the first hand, the real axis is divided in the following intervals 

+ +.../ +/ +...+ / (B.3) 

1 ,u) J,.u) L .L.u) 



Vi J Vk 



"3 



(in strip 2 there is no need to consider intervals with £). The integration variable can be changed 
to be t = tj, with dt = t'-dx and in this new variable it is apparent that all the contributions 

corresponding to the intervals (|B.3|) vanish because tj{ytp) = 0, except the first and the last. This 
leads to 

S] =i M+ ™\tM i+t \-'^m. (b.4) 



12 



Along the flow (£ EM), the required endpoints of the integral take the following values: 



1 + 

tj(+oo) = 2 cos A = — - — , ^ R5 j 
ti(-oo) = si log gx (-00, £) = -1, ^(-00) = 0. 

The integral itself is recognised to be a sum of Rogers dilogaritms that, at the previous endpoints, 
can be exactly computed. Actually only the sum is required: 

Si + S 2 = ^vr 2 . (B.6) 
15 

On the other hand, one can substitute in (|B.2|) the expressions for €j, e'j given by the TBA equations 
(|4,3( I4.4jl . In doing this, all terms involving the convolution have cancelled because of the symmetry 
of the kernel, K(x) = K(—x). Among the remaining terms, those involving the derivatives V- are 
integrated by parts leading to 



S 1 + S 2 



dx £ 1 + £ * - 4e-L 2 +2 / d X f Ll . (B.7) 

k=1 smh(y K k - x) k=1 smh(y k - x) 91 



The first term on the right hand side enters the expression for the energy (|B.1|) . Using the result 
(jB.6|) one obtain, for all £ E K: 

1 1 T/o 2 m i • m 2 r+00 j -/ 

B (f) = -iM + — i— '■EE J f , -sFi; ,og< "^ , -2'0 + y«, 4?| L " (R8) 

i=l fc=l fe=l J ~°° y 

This is the espression for the energy with the most explicit dependence on the boundary term g\ 

so that is well suited to the purpose of taking the limit £ — > +00. Using the explicit expression for 

(2) (2) 
gi and remembering that, along the flow, y^ = —00 if and only if = n 2 , one obtains 

i m2 7r i m2 y^ + £ 7r 1 

^-^E^rf ) - = E lQ g tanh (^— - i£) = -2^,n 2 - ( B - 9 ) 

fe=l k=l 

This term is the new feature of the flow under examination with respect to the cases with single 
character. In spite of the fact that lim^^ +00 g[ = 0, the term under the integral also gives a 
contribution because of the divergence L\{— 00) = logO = —00 (|B.5|I . From the TBA equations the 
following asymptotic behaviour is obtained, for all fixed values of £ 

ei(x) ~ -logtanh^-L (B.10) 

To get this from (|4,3[) . one considers that y^ 3> —00 because of the infinite number of 2-strings 
appearing in the scaling limit p. In addition to that, L 2 (— 00) = was used. Actually, the 
following term will be used: 

, _ Sl e-^e[ ^ i( 1 + tanh £±l). ( B . n) 



This asymptotic behaviour in the variable x is relevant in the integral term in (|B.8|) . when £ — ► +00. 
This is apparent from the following explicit expression: 

„ r+°° dx g\ T f +oc it, [Tl f + °° dr, , i + { lr , ^ . 
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where an integration by parts has been performed and the corresponding boundary terms vanish. 
One has: 



1 --/ —logltanh— Kl + tanh — 



oo 

= — / — r log I tanhxlfl + tanhx) 
J-oc 47r 

/■ + °° <£c , , f +c ° dx , 

= — / — 7 log tanh a; = — / — ^logtanhx (B.13) 

./-oo 47r JO 27r 

where on the first line the asymptotic behaviour was used, in the second line a redefinition of x 
shows that the integral is independent on £, on the third line a term vanish because it is odd under 
the exchange x — > — x; the remaining integrand is even under the same exchange and this leads to 
the last equality. With a new change of variables, this definite integral is transformed in 

f 1 dt 1 , 1 - 1 1 , 
T=- I — ~ — log = — . B.14 

J 4vr 2 t 1 + 1 16 y J 

The actual value can be explicitly computed or obtained from standard mathematical tables of 
integrals; it admits also an expression in terms of Rogers dilogarithms. In conclusion, from (|B~ 
IB .91 IB.14|) one has the energy at the IR critical point (1,1) © (3, 1) 



„ 7 m 1 Cm 1 . v^v^t(«) 

E = -™\ + 4 " o 6 l?\n 2 + E E 4 } (B.15) 



i=l k=l 



that is precisely the expression (|2.7|) given in Section [2 The first term is the central charge 



contribution ~£z- 



As a final comment, the expression ()B.8|) can be generalized to all the flows with the obvious 
addition of terms depending upon gi and a possibly non-zero term that takes into account the 
subtleties related to t±(— oo). 



Appendix C. Rule for the integration constants sj. 

The braid limit used in the first line of (|B.5|) can be used to obtain the general selection rule for 
the variables Sj. They appear in the derivation of the TBA equations as integration constants (see 
for details). Then, taking the limit x — > +oo in the equations themselves, one can use the braid 
limit to express the asymptotic values of ej, Lj. This leads to: 

8j = {-l) m i , f G R, £ = +oo = IR. (C.l) 

An implicit assumption was used, that all the 1-strings are in a finite position, < +oo, condition 
that is always satisfied except at the UV point. If this is violated, the actual value of the braid 
limit can change and the previous argument doesn't hold. The correct UV rule can be obtained by 
the previous one observing that if along the flow one has c/j ^ 1, the content of zeros changes in 
the strip 3 — j so that: 

Sj = (-iy^+\ £ = _oo = UV. (C.2) 

The selection rules proven here hold for all the flows and all the conformal boundary conditions. 
Of course, they contain the cases given in and, in this sense, this Appendix completes the 
exposition given there. 
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